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Abstract 

In  this  paper  several  previously  developed  theories  for  predicting  liftoff  heights  of 
turbulent -jet  flames  are  examined  and  compared.  The  five  theories  are:  a  turbulent  flame 
speed  theory,  a  cascade  mixing  time  theory,  a  strain-rate  theory,  a  scalar  dissipation-rate 
theory  and  an  empirical  correlation.  Each  theory  uses  far-field  scaling  of  isothermal, 
turbulent  jets.  Through  a  systematic  algebraic  development,  the  predictions  of  the 
theories  are  presented  in  a  common  form.  They  are  then  compared  for  their  ability 
to  predict  the  experimentally  observed  linear  dependence  of  liftoff  height  on  nozzle  exit 
velocity.  Following  this,  their  ability  *o  correctly  predict  the  relative  liftoff  heights  of 
different  fuels  is  examined.  Finally,  the  theories  are  extended  to  predict  blowout  and 
their  predictions  are  compared  to  experimental  results  of  several  fuels. 


List  of  Symbols 
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a  velocity  gradient  at  extinction  in  a  counterflow  diffusion  flame 
c,  constants,  i  =  1,  2,  3,  ... 

cv  constant  with  dimensions  of  kinematic  viscosity 

Cp  specific  heat  at  constant  pressure 

d  nozzle  diameter 

d*  momentum  diameter,  d(pQ/Poo) 

D  molecular  diffusivity 

ft{rj)  function  of  r/  as  defined  by  a  particular  equation,  i  =  1,2,3,  ... 

Fn(r ])  ‘‘correction  functions”  to  make  theories  match  experimental  results,  n  =  t.  m,  s,  d  and  e. 

h  liftoff  height 

k  thermal  conductivity 

l  integral  length  scale 

r  radial  coordinate 

5  maximum  laminar  flame  speed 

u  mean  velocity 

u'  rms  velocity 

U0  nozzle  exit  velocity 

x  axial  distance 

V  fuel  mass  fraction 

Y0  initial  fuel  mass  fraction  (balance  is  air) 

T]  r/x 

k  thermal  diffusivity,  k/(pCp) 

v  kinematic  viscosity 

p  density 

(7  strain  rate 
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tc  fuel  chemical  time 

Tf  cascade  mixing  time 

X  scalar  dissipation  rate,  2£>|VFj2 

Other  symbols 
(.)0  nozzle  exit  conditions 

(,)oo  reservoir  conditions 

(,)3t  stoichiometric  conditions 

(.)/  conditions  at  maximum  laminar  flame  speed 
(-)stab  assumed  stabilization  conditions,  either  (,)st  or  (.)/ 
{■)qu  quenching 

(.)d  centerline 

blowout 

(.)t  turbulent  flame  speed  theory 

(.)m  cascade  mixing  time  theory 

(.)3  strain-rate  theory 

(,)d  scalar  dissipation-rate  theory 

( .  )e  empirical  correlation  theory 

(.)„  name,  representing  t ,  m,  s,  d  or  e 

(.)  time  average 
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1.  Introduction 

Several  theories  using  far-field  scaling  of  isothermal,  turbulent  jets  have  been 
developed  to  explain  the  dependence  of  the  mean  liftoff  height  above  the  nozzle  exit, 
/?,  on  the  flow  parameters.  By  appropriately  grouping  the  parameters  in  each  theory,  the 
dependence  of  h  on  the  nozzle  velocity  U0  momentum  diameter  d*  initial  mass  fraction 
of  fuel  Y0  and  chemical  and  transport  properties  of  the  fuels  is  brought  out. 

The  theories  examined  here  are:  the  turbulent  flame  speed  theory  of  Kalghatgi  [1]. 
the  cascade  mixing  time  theory  of  Broadwell  et  al.  [2],  the  strain-rate  theory  of  Miake- 
Lye  and  Hammer  [3],  the  scalar  dissipation- rate  theory  of  Peters  and  Williams  [4]  and  a 
theory  developed  by  Pitts  [5].  Each  of  the  above  theories  is  treated  in  a  similar  manner. 
First,  through  algebraic  manipulation,  an  equation  is  derived  involving  h,  U0 ,  Y0,  fuel 
parameters,  similarity  functions  of  rj  =  r/x  for  velocity,  concentration,  etc.  and  what 
will  be  called  a  correction  function,  Fn(r}),  which  is  needed  to  make  predictions  match 
experimental  results.  (Ideally,  Fn(r/)  =  1  for  each  theory.)  Next,  the  necessary  shape 
of  Fn(r})  to  make  the  theory  match  results  is  found.  Finally,  setting  Fn(r])  =  1,  r/  is 
eliminated  from  the  equation  and  the  true  predictions  of  the  theory  are  given. 

The  remainder  of  the  paper  is  organized  as  follows.  First,  experimental  results 
of  lifted  flames  are  summarized  and  assumed  scaling  laws  of  isothermal  jets  are  given. 
Next,  each  of  the  five  theories  is  analyzed  as  described  in  the  preceding  paragraph.  In 
the  Comparisons  section  the  similarities  and  differences  of  the  theories  are  summarized. 

2.  Experimental  results  and  assumed  similarity  profiles 
2.1  Results  for  lifted  jet  flames 

Kalghatgi  [1]  has  found  that  for  pure  fuel  issuing  from  a  nozzle  h  is  proportional  to 
U 0  and  independent  of  d  for  a  given  fuel,  once  h  reaches  ~  10d  In  non-dimensional  terms 


this  can  be  written  as 


h_ 

A 


(1) 
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=  T„ 


,Uo 


where  t'c  is  a  chemical  time  dependent  of  the  fuel.  For  small  dilutions  (F0 
h/U0  oc  ( 1  / }  o)  [3]  and  so  Eq.  (1)  can  be  written  as 


1) 


A  =  A2A 

d*  ~ Tc  d *  Y0  ' 


(2) 


For  larger  dilutions  h/U0  was  not  found  to  be  proportional  to  1  /Ya,  although  this  may 
be  due  to  near-field  effects. 

This  liftoff  behavior  persists  until  the  flame  blows  out  at  a  maximum  velocity  Ui, 
which  is  proportional  to  d,  with  the  constant  of  proportionality  depending  on  the  fuel. 


2.2  Assumed  scaling  laws  of  isothermal  turbulent  jets 

The  self-similar  jet  profiles  which  are  used  in  the  theories  are  the  following: 


u  =  c\U0— f\(r]) 

X 

(3) 

h{-n)  =  e-c'y 

(3.1) 

Y  =  c2Y0 — f2(r]) 

X 

(4) 

fi(v)  =  ^ 

(4.1) 

U'  =  C3U0  —  f3{7]) 

X 

(5) 

/3(d)  =  t~c^ 

(5.1) 

d* 

0  =  c4Li0~f4(j] ) 
xl 

(6) 

U  (Tl)  =  l~cW 

(6.1) 

=  c5  — )4/o(^) 

a  x 

(7) 

/  =  CGX 

(S) 
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The  above  Gaussian  functions  (Eqs.  (3.1),  (4.1),  (5.1)  and  (6.1))  are  used  only  for 
estimating  true  predictions  of  the  theories.  Functions  of  any  shape  may  be  used  in  place  of 
the  Gaussian  functions  at  the  expense  of  greater  complication  in  stating  the  final  results. 

3.  Brief  development  of  each  of  the  five  theories 

In  this  section  the  five  theories  (turbulent  flame  speed,  cascade  mixing  time,  strain 
rate,  scalar  dissipation  rate  and  the  empirical  theory)  are  examined.  The  descriptions  do 
not  exactly  follow  the  authors’  developments  for  each  of  the  theories.  This  is  done  in 
order  to  bring  out  the  similarities  and  differences  in  the  theories.  Major  deviations  from 
the  authors’  developments  are  noted. 

3.1  Turbulent  flame  speed  theory 

The  turbulent  flame  speed  theory  assumes  that  the  mean  velocity  u  equals  the 
turbulent  flame  speed  ut  at  the  liftoff  location  specified  by  both  axial  x  and  radial  r 
coordinates.  Also,  the  liftoff  location  is  assumed  to  lie  along  the  contour  of  maximum 
laminar  burning  velocity,  Y  =  Yi.  Ut  is  assumed  to  depend  on  the  flow  properties  as 

-T  =  c'7(u'l/ui)1/2  =  c7(u'x/ui)1/2.  (9) 

There  is  some  experimental  support  for  this  expression,  although  considerable  scatter 
exists  [6]. 

Kalghatgi  [1]  uses  dimensional  analysis  in  his  development  and  eventually  arrives 
at  an  expression  for  thp  turbulent  flame  speed.  The  current  description  starts  with  the 
assumption  of  the  turbulent  flame  speed  relation,  Eq.  (9),  and  arrives  at  an  equation  for 
h.  Its  development  is  different  from  that  of  Kalghatgi  in  this  way  and  in  the  use  of  the 
correction  function  described  below. 

Rather  than  the  rigid  requirement  u  =  ut  at  the  liftoff  position,  in  the  present 
discussion  we  admit  that  jet  turbulence  is  not  homogeneous  and  isotropic  and  that  there 
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will  be  differences  at  different  radial  locations  rj.  The  correction  function  for  this  theory 
is  denoted  as  Ft(rj)  and  we  require 


u  =  Ft(rf)ut , 


(10) 


where  Ft(r/)  =  1  for  the  turbulent  flame  speed  theory  to  be  exact.  Squaring  Eq.  (10). 
using  Eq.  (9),  and  setting  x  =  h  yields 


«2  = 

Vl 


Using  Eqs.  (3)  and  (5), 


ci  fi(v) 
cs  fsiv) 


u'  = 

Vl 


Using  Eq.  (4)  and  the  assumption  Y  =  Yi  to  eliminate  d*/h, 


c2c3  Yo  /2(r/)/3(r?) 


vi 


Finally,  solving  for  h , 


i  Yl  u‘  tt  /i2(7?)/(/2(^)/3(^)) 

n  ~  ct,\TF  —  Uo - 


Y0  S2 


F?(V) 


(ID 


where  the  constants  have  been  absorbed  into  c^j. 

Some  important  points  to  note  from  Eq.  (11)  are  that  at  a  fixed  rj,h~  1  /Y0,  h  ~  U0. 
and  h  does  not  depend  on  d* .  These  relations  will  turn  out  to  be  true  for  all  of  the 
theories  in  this  paper. 

There  are  two  remaining  issues  to  address  regarding  the  turbulent  flame  speed 
theory.  The  first  is  whether  /<(r?)  =  ~  constant  as  required  for  agreement  with 

experimental  liftoff  results  that  there  is  a  linear  relation  between  h  and  Ua  for  a  given  fuel. 
By  determining  the  actual  shape  of  /<(r/)  based  on  experimental  measurements  of  fi(rj), 
.^2 ( r/ )  and  /:}(r?).  the  required  correction  function  Ft(rj)  can  be  determined.  Alternatively. 
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once  ft  ( Jj )  is  found.  Ft(rj)  can  be  set  to  1  to  show  the  true  turbulent  flame  speed  theory 
prediction  without  correction. 

The  second  issue  is  whether  the  theory  correctly  predicts  the  relative  liftoff  heights 
of  different  fuels.  This  is  measured  by  whether  the  factor  Yiui/S2  in  Eq.  (11)  agrees 
with  experimental  results  of  slopes  of  h  versus  UQ  lines  for  various  fuels  (r'  in  Eq.  (1)). 
Addressing  this  second  issue  is  postponed  until  the  Comparisons  section  where  all  five 
theories  are  compared. 

Returning  to  the  first  issue,  experimental  data  is  used  to  approximate  the  functions. 
/o( !])  and  /;d?/)  with  Gaussian  profiles  as  indicated  in  Eqs.  (3.1).  (4.1)  and  (5.1). 
For  mean  velocity.  c\  =  93.7  from  data  of  Chen  and  Rodi  [7].  For  fluctuating  velocity. 
r't  =  40.5  is  found  from  curve  fitting  data  of  Wygnanski  and  Fiedler  [S]  at  Y /Y,-i  =  1/2. 
This  curve  fit  turns  out  to  be  about  30c/c  high  near  the  jet  edge.  ??  =  0.2.  but  will 
nevertheless  be  used  here  for  simplicity.  Finally,  for  concentration,  c'2  =  57  [7.  9].  This 
results  in 


ft(n)  =  e-(2c;-c'2-ci),1  =  e-cW  =  t-wn\  (12) 

The  ratio  of  /,  at  rj  =  0.1  and  0.2  (the  approximate  range  of  values  encountered 
experimentally  in  lifted  flames)  is  e-4o('2  *  =  0.2G.  Therefor.  Ft(i])  would  need  to  be 

%  (  _  and  the  ratio  u/vt  would  need  to  vary  by  a  factor  of  ~  4  to  make  the  turbulent 
flame  speed  theory  match  experimental  results. 

If  Ft( )/)  is  set  identically  equal  to  1.  the  predictions  of  a  true  turbulent  flame  speed 
theory  can  be  found.  Multiplying  Eq.  (11)  by  Y)/(Y0d*), 

iiA  _  .  a 

Yn  cl*  C,A  Y0  ]  S2  c l*  -  ’  11  ' 


Solving  for  //  from  Eq.  (4) 


\  i  2  v  i  L  „  _  j 

C,AiY0)  S2  d-ftif'2  ( 


1AA 

c-2  F0  d* 


))- 


(13) 
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Using  the  Gaussian  functions  results  in  the  power  law  dependence. 


AlA  =  c 

y; d*  tM[Y0}  s2  d *  ’ 


Yl  .  2  Vl  Uo 


(14) 


The  exponent  Et  =  _  )-j-  =  —1.73.  Because  Et  <  0,  the  turbulent  flame  speed  theory 

1  Ct/C2 

actually  predicts  a  decreasing  liftoff  height  with  increasing  nozzle  velocity  for  a  fixed  d* 
and  fuel. 


3.2.  Cascade  mixing  time  theory 

In  a  paper  devoted  primarily  to  blowout  of  flames,  Broadwell  et  al.  [2]  proposed  that 
lifted  flames  will  be  stabilized  at  a  distance  such  that  the  mixing  time,  which  is  controlled 
by  a  turbulence  cascade  time,  tj  ~  x/uci,  is  equal  to  a  chemical  time  for  the  fuel.  rcm . 
Their  theory  is  one  dimensional  because  tj  is  assumed  to  have  to  no  dependence  on  r/. 
rc  m  is  assumed  to  be  proportional  to  k/S 2  in  accordance  with  simple  one-dimensional 
laminar  flame  speed  theories  (e.g.,  Glassman  [10]).  A  discussion  of  what  value  of  k  should 
be  used  is  postponed  until  the  Comparisons  section. 

For  the  present  description,  a  correction  function  Fm(Tj)  will  be  added  to  the  theory 
as  was  done  with  the  turbulent  flame  speed  theory.  Using  Eq.  (3)  the  liftoff  requirement 
becomes 

=  C,„  Jj.  (15) 

Assuming  the  flame  stabilizes  along  the  mean  stoichiometric  contour  (an  assumption 
that  was  not  necessary  in  the  original  one-dimensional  formulation)  and  using  Eq.  (4), 


h  = 


Xli  —  TT  1/h(rl) 
CmJ  y;  s2  °  Fm ( rj )  ■ 


(16) 


As  with  the  turbulent  flame  speed  theory,  at  constant  i),  h  ~  Un ,  h  ~  1/U„  and  h  is 
independent  of  d* . 


To  find  the  true  predicted  liftoff  height,  Fm(r} )  is  set  to  1  in  Eq.  (16)  to  get 


Yst  h _  /  Yst  ,2  K  -t  /  r  r  \ 

Y  d*  ~~  V  C2  rl*  ! fevl)- 


S2  d * 


(17) 


Noting  the  similarity  to  Eq.  (13),  we  can  write 


Yst_h_ 

Y0  d* 


,2  K  \  1  /2 


'  Cma((  )  ^  d 


(IS) 


The  Gaussian  shape  of  f2(r) )  was  not  used  in  deriving  Eq.  (IS)  because  no  other  functions 
were  involved,  and,  the  exponent,  Em  =  1/2,  is  exact.  In  fact,  Eq.  (18)  can  easily  be 
derived  from  the  strictly  one-dimensional  theory. 


3.3  Strain-rate  theory 

A  theory  very  similar  to  the  above  mixing  time  theory  is  the  maximum  strain-rate 
theory  of  Miake-Lye  and  Hammer  [3].  Here,  however,  the  liftoff  condition  is  that  the  large 
scale  strain  rate  a  cannot  exceed  a  maximum  strain  rate  crcrit  ~  S2/k.  Also,  a  is  allowed 
to  have  a  radial  dependence,  hence  making  this  a  two-dimensional  theory.  The  flame  is 
assumed  to  lie  along  a  contour  on  which  stoichiometric  conditions  are  most  probable.  To 
approximate  this  condition,  it  is  assumed  that  Y  =  Yst- 

Using  a  correction  function  Fs(r 7)  the  liftoff  requirement  becomes 


d*  S2 

a  =  c4U0  — /4(r ))  —  cn—Fs(T]). 


x 


K 


Solving  for  h, 

Yst  K  rr  U(v)/f2(ri) 

which  is  very  similar  to  Eq.  (16).  The  only  difference  is  that  the  radial  strain  function 
/4(r/)  appears.  Therefore,  to  match  experimental  liftoff  results,  Fg(r])  =  fi(r/)/ fair)). 

Unfortunately,  U(t])  has  not  been  experimentally  measured.  Nor  was  it  completely 
defined  in  Ref.  [3]  where  the  authors  described  it  as  the  large-scale  strain  rate  between 


large  vortical  structures.  One  possibility  is  to  use  ^  for  a.  This,  however,  will  not  give  a 
Gaussian  profile.  Instead,  in  the  present  paper,  a  Gaussian  profile  is  used  which  matches 
at  r)  =  0.1  and  77  =  0.2  to  get 


c\  =  c\  -  log(. 2/.l)/(.22  -  .l2)  =  71. 


This  makes  fs(rj)  =  fi{r])lf2(r))  =  e  ^  where  c's  =  c\ 


c'2  =  14. 


The  true  predictions  of  the  theory  are  then 

•) 


hiJL-r  /(i— c; / c'2) 

Y0  d *  ~  S’2U 


Y0  '  S2  d' 


-)1 


(20) 


The  exponent  E3  =  1/(1  —  c's/c'2)  =  1.3,  although,  due  to  the  approximation  in 
determining  this  number  should  be  treated  with  some  skepticism. 


3.4  Scalar  dissipation-rate  theory 

Peters  and  Williams  [4]  have  proposed  that  the  liftoff  height  is  determined  by 
extinction  of  thin  flamelets.  Through  an  asymptotic  analysis  they  find  that  the  scalar 
dissipation  rate  y  is  the  parameter  which  determines  whether  a  flamelet  will  be  quenched. 
The  criterion  for  flame  stabilization  is  that  the  probability  of  \  <  Xqu  at  stoichiometric 
conditions  is  greater  than  some  minimum  probability,  where  \qu  is  a  limiting  scalar 
dissipation  rate  for  the  particular  fuel.  Using  percolation  theory,  this  requirement  is 
simplified  to  \st  <  Xqu ■  In  the  present  formulation,  we  will  again  add  a  correction 
function  Fd(rj)  and  require 

Xst  =  XquFdin)-  (21) 

Here,  \3t  is  the  average  scalar  dissipation  rate  conditional  on  Y  =  Yst.  In  order  to 
make  use  of  experimental  measurements  in  jets,  must  be  related  to  the  unconditional 
average  dissipation  rate  F°r  this  we  assume  that  Xst  =  xfi^st/Yci,^)^  and  that 
Y  =  Y3t.  Then  \st  =  xf7(f2(v),T])-  Combining  this  with  Eq.  (21), 

V  =  Xstlh(ri)  =  Xqu  Fd(r])  /  f~(r}). 


Now,  assuming  /V (77)  is  not  a  strong  function  of  77  (or  equivalently  absorbing  it  into  Fd(T 7)) 
we  arrive  at  the  condition 


X  =  c,2XquFd’l) 


(22) 


for  determining  the  liftoff  height. 

After  that  somewhat  complicated  introduction,  we  are  no  ready  to  proceed  using 
Eqs.  (22),  (7)  and  (4)  as  with  the  previous  theories  to  get 
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''YoX^9  Fd(v) 

Like  the  previous  theories  at  constant  77,  h  ~  U0,  h  ~  1  jY0  and  h  is  independent  of  d* . 
Unlike  the  previous  theories  an  additional  Y2t  appears  in  the  result,  but  Y2t/xqu  is  simply 
a  chemical  time  relating  the  relative  liftoff  heights  of  different  fuels. 

For  estimating  the  shape  of  /s(^),  results  of  Dowling’s  [9]  point  measurements  of 
concentration  are  used.  Using  Taylor’s  hypothesis  along  with  the  mean  off-centerline 
velocity,  he  states  that  “the  mean  value  of  [the  estimated  scalar  dissipation  rate]  is  almost 
independent  of  [77]  ...”.  Therefore,  ^5(77)  is  approximated  as  a  constant,  and  Fd{r 7)  would 
need  to  be  equal  to  fd(r])  =  /s(?7)//| (77)  ~  e-171??2  for  the  scalar  dissipation-rate  theory 
to  match  experimental  results.  Peters  and  Williams  approximated  \  as  U0d*(^)2  which 
makes  /"s (77 )  ~  r)2e~2c'iin  if  the  Gaussian  profile  for  /2( 77)  given  by  Eq.  (4.1)  is  used.  The 
current  discussion  will  assume  that  ^(77)  sa  constant ,  however. 

Setting  Fd{rj)  =  1  in  Eq.  (23)  results  in 


Y0  d 


Xqu 


(24) 


where  c'd  =  —  3c'2  and  therefore  Ed  =  1/(1  —  c'd/c2)  =  1/4. 

It  should  be  pointed  out  that  Peters  and  Williams  have  not  stated  that  their  theory 
should  be  applicable  to  diluted  fuels  where  Y0  <  1.  The  idea  that  the  quenching 
scalar  dissipation  rate  \qu  is  a  constant  is  based  in  part  on  the  assumption  of  boundary 


conditions  of  the  flamelets  being  pure  fuel  on  one  side  and  pure  oxidizer  on  the  other. 
However,  for  comparison  with  other  theories  in  this  paper,  the  effect  of  Y0  has  been 
included  with  the  assumption  that  Xqu  would  not  change. 

One  other  point  to  be  made  concerns  the  constant  C5  in  Eq.  (7).  Although  Eq.  (7) 
is  expected  to  be  correct  in  the  limit  Re  — >  00,  at  Reynolds  numbers  of  5000  and  16000, 
Dowling  finds  C5  =  54  and  32  respectively.  This  would  cause  even  further  disagreement 
between  experimental  results  and  the  scalar  dissipation-rate  theory  in  this  range  of 
Reynolds  numbers. 

3.5  Empirical  correlation  proposed  by  Pitts 

Pitts  [5]  has  proposed  a  correlation  for  predicting  liftoff  heights  which  is  based  on 
jet  scaling  (with  some  caveats).  He  proposes  that  flames  stabilize  along  the  contour  of 
maximum  laminar  burning  velocity  (Y  =  F/)  at  a  location  where 

a  =  (25) 

Ci/  d 

The  first  difficulty  with  the  correlation  is  that  c„  is  a  constant  with  dimensions  of 
kinematic  viscosity,  independent  of  fuel  and  oxidizer,  which  cannot  physically  be  correct. 
A  possible  fix  to  this  would  be  to  replace  cv  with  tq,  although  this  would  change 
predictions  slightly  for  different  fuels. 

The  second  difficulty  concerns  the  factor  Y2  in  Eq.  (25).  Whether  this  should  be 
replaced  by  ( Y[/Y0 )2  or  ( Yi/Y0)Yi  or  left  as  Y2  when  Ya  <  1  is  not  stated  by  Pitts. 
The  correlation  was  developed  to  match  data  for  pure  fuels  only,  and  so  this  was  not 
a  consideration.  In  the  present  discussion  Y2  will  be  replaced  by  (Yi/Y0)Y)  because 
the  other  two  choices  result  in  ustab  depending  on  x/d or,  through  simple  algebraic 
manipulation,  on  Y0.  It  is  assumed  in  jet  far-field  scaling  that  the  flow  is  independent  of 
these  parameters. 
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Replacing  Y'2  with  Y )(Y)/Y’0)  in  Eq.  (25)  and  adding  a  correction  function  Fe(rj). 

u  =  -^xS2Y,(Yl/Y0)Fe(r]).  (26) 

c„  d * 

Using  Eq.  (4)  to  eliminate  (x/d*)(Yi/Y'0)  better  illustrates  the  assumed  dependence  of  u 
on  the  flow  variables: 

u/S  =  ciYixS/c^MriF'it]).  (26.1) 


It  is  difficult  to  justify  theoretically  this  relationship,  especially  with  the  factor  Y).  but 
at  least  there  is  no  remaining  x/d*  which  would  put  it  out  of  the  realm  of  jet  scaling 
theories. 

Using  Eq.  (3)  in  Eq.  (26), 
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where  Yi  has  been  explicitly  left  in  the  numerator  and  denominator.  For  a  linear 
relation  between  UQ  and  ft,  Fe(r/)  must  equal  fe{r])  =  /i(t7)//| (v)  =  e~c‘v7 ,  where 
c'e  =  c\  —  2c'2  =  —20.  Pitts  effectively  used  such  a  correction  function  by  using  a  slightly 
higher  value  of  c\  than  had  been  found  in  turbulent  jet  experiments. 

The  true  predictions  of  the  theory  are 


Yh_ 

Y0  d* 


=  c 

In 


C"  )E, 

Y)S 2  d*  ’  ’ 


(28) 


where  Ee  —  1/(1  —  c'e/c2)  =  0.74. 


4.  Comparison  of  theories 

Each  of  the  theories  investigated  produces  a  prediction  of  the  form 

h  =  cn,l(Ystab/Y'0)Tc,nU0fn(i1 ),  (29) 

where  Y\tab  is  the  assumed  stabilization  concentration,  Yrst  for  the  mixing  time,  strain-rate 
and  scalar  dissipation-rate  theories  and  Y\  for  the  turbulent  flame  speed  and  empirical 


theories.  rCi„  is  a  fuel-dependent  chemical  time  and  fn (rj)  is  a  function  which  depends  on 
jet  similarity  functions.  The  above  relation  is  a  direct  result  of  the  fact  that  the  proposed 
stabilization  conditions  of  each  of  the  theories  can  be  expressed  as 

- - 7-c,n/n.2(7?), 

Ucl 

Y  =  Ystab , 


where  /n> 2(r?)  =  fn{v)f2(T])-  By  eliminating  r\  from  Eq.  (29)  using  Eq.  (4),  the  liftoff 
prediction  can  be  stated  as 


}  stab  h 


/n,3(( 


(30.1) 


Approximating  /„(r;)  and  ^(r?)  with  Gaussian  functions  results  in  a  simple  power  law 
form  for  the  function  /n)3 


^  stab  h 

~Yo~~d* 


(30.2) 


This  section  is  organized  into  four  subsections.  First,  predictions  for  a  fixed  fuel  (and 
fixed  chemical  and  transport  parameters)  are  compared  by  examining  /n(^)-  Second,  the 
necessary  parameters  for  comparing  different  fuels  via  different  expressions  for  Ystabrc,n 
are  tabulated.  The  third  subsection  uses  these  parameters  to  compare  experimental 
results  for  liftoff  heights  with  predictions  for  different  fuels.  The  last  subsection  addresses 
blowout  predictions. 


4.1  Predictions  for  a  single  fuel 

Table  1  lists  the  following  properties  of  each  theory:  expression  for  /n(r/)  in  terms 
of  similarity  functions  /j,  /2,  etc.;  the  approximate  value  of  the  c'n  when  fn(ri)  is 

/  2 

approximated  as  e~CnV  ;  the  the  exponent  in  Eq.  (30.2),  En  =  1/(1  —  c'n/c'2)-,  and  the 

* 


expression  for  YstabTc 


Of  the  above  theories,  only  the  turbulent  flame  speed  theory  predicts  a  decreasing  h 
with  increasing  U0 ,  as  evidenced  by  En  <  0.  The  scalar  dissipation-rate  theory  has  the 
largest,  negative  c'„  and  consequently  predicts  a  dependence  of  l7®  25.  Other  theories  have 
a  more  nearly  linear  dependence  in  agreement  with  experimental  results.  Due  to  lack  of 
data  on  the  strain-rate  profile,  the  results  for  the  strain-rate  theory  are  less  certain  than 
those  of  the  other  theories. 


Table  1 

Summary  of  Theory  Predictions 


Theory 

fniv) 

< 

En 

h  stable, n 

turbulent  flame  speed 

/iV(/2/3) 

90 

-1.73 

Y^vi/S2) 

cascade  mixing  time 

l//2 

-57 

1/2 

Yst(K/S2) 

strain  rate 

U/h 

14 

1.3 

Yst(*/S2) 

scalar  dissipation 

hlfl 

-171 

0.25 

Yst(Ys2t/Xq  u) 

empirical  correlation 

f\!fl 

-20 

0.74 

Yi(cJ(YiS2)) 

4.2  Properties  for  evaluation  of  chemical  times 

In  Table  2,  values  for  K,<,  Y)  and  5  are  taken  from  Kalghatgi  [1].  In  accordance  with 
turbulent  flame  speed  correlations,  i//  is  evaluated  at  the  unbumed  condition,  assumed 
to  be  300  I\. 

Evaluation  of  k  is  somewhat  more  difficult.  Broadwell  et  al.  used  for  a  value  for 
k  that  of  air  at  2000/i  for  all  the  fuels  because  the  hydrocarbons  have  similar  flame 
temperatures.  (Hydrogen’s  temperature  is  somewhat  higher.)  Miake-Lye  and  Hammer 
used  k  of  stoichiometric  unbumed  mixtures  at  2200  K,  although  this  required  considerable 
extrapolation  for  the  thermal  conductivities  of  some  of  the  components.  More  detailed 
inspection  of  the  laminar  flame  speed  theories  reveals  that  p  should  be  evaluated  at 
the  cold,  unburned  condition.  The  theory  of  Mallard  and  Le  Chatelier  [10]  requires 
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Cp  to  be  evaluated  at  a  temperature  between  the  flame  temperature  and  the  unbumed 
temperature,  and  k  to  be  evaluated  near  the  flame  temperature  for  the  most  physical 
interpretation  of  the  theory.  The  more  complex  theory  of  Zeldovich,  Frank-Kamenetskii 
and  Semenov  [10]  requires  that  Cp  and  k  be  constant.  Therefore,  in  the  Table  2,  p  and  C'p 
are  evaluated  at  stoichiometric,  unburned  conditions,  and  k  is  evaluated  at  the  adiabatic 
flame  temperature.  Because  reliable  data  for  k  of  the  products,  H20  and  CO2,  could  not 
be  found,  k  for  N2  at  the  flame  temperature  is  used. 

Values  of  \ q u  are  obtained  from  the  equation  of  Peters  and  Williams 

Xqu  =  Xst  «  ^aY^[erfc-1(2Yat)}2,  (31) 

where  the  velocity  gradient  at  extinction  in  a  counterflow  diffusion  flame  a  is  taken  from 
the  data  of  Tsuji  and  Yamaoka  [11]  for  CH4  and  CgHg  and  from  the  computations  of 
Dixon-Lewis  et  al.  [12]  for  H2.  No  value  has  been  found  for  C2H4. 

Also  included  in  Table  2  are  experimental  results  for  the  slopes  of  the  h  versus  Ua 
lines  from  the  data  of  Kalghatgi  [1].  The  final  column  gives  data  on  blowout  in  the  form 
of  d* /Ub  taken  from  Broadwell  et  al.  who  used  the  data  of  Kalghatgi  [13].  These  last  two 
columns  will  be  compared  with  predictions  in  sections  4.3  and  4.4. 

Table  2 

Fuel  properties  and  experimental  results  for  liftoff  and  blowout 


Fuel 

Y,t 

Yi 

5 

m/s 

v\ 

10  ~6m2/s 

K  St 

10~6m2/s 

Xqu 

■s-1 

h/U0 

ms 

d*  /Ub 

ps 

ch4 

0.0549 

0.058 

0.39 

15.9 

106 

4.9 

2.42 

46.6 

c2h4 

0.0632 

0.073 

0.75 

15.1 

106 

— 

0.672 

19.4 

C3Hs 

0.060 

0.068 

0.45 

14.6 

102 

6.1 

1.91 

51.2 

h2 

0.0284 

0.076 

3.06 

30.1 

113 

76.1 

0.030 

0.178 
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4.3  Comparison  of  predicted  and  measured  liftoff  heights  of  different  fuels 

Several  points  can  be  noted  about  the  expressions  for  Ystabrc,n  given  in  Table  1. 
All  contain  a  1/5*  term,  except  the  scalar  dissipation-rate  theory.  The  expressions  for 
the  mixing  time  and  strain-rate  theories  are  identical.  Also  the  empirical  correlation’s 
expression  would  agree  with  that  of  the  turbulent  flame  speed  theory  if  cvjY\  were 
replaced  with  v\  in  Eq.  (2G.1). 

Using  the  properties  in  Table  2,  the  relative  predicted  liftoff  slopes  (YstabTc,n)  of  the 
fuels  are  evaluated  and  tabulated  in  Table  3.  For  ease  of  comparison,  the  values  have 
been  normalized  by  the  value  for  methane,  as  have  the  experimental  results  for  h/U0.  In 
addition,  the  deviation  of  the  prediction  from  the  experimental  result  is  listed  in  percent. 
All  of  the  theories  correctly  predict  the  relative  liftoff  heights  of  different  fuels  to  within 
±  ~  30%  except  for  the  turbulent  flame  speed’s  prediction  for  hydrogen. 

Table  3 

Comparison  of  measured  and  predicted  relative  liftoff  heights  for  several  fuels 


Fuel 

h/U0 

turb.  f.s. 

mix.  time 

strain  rate 

scalar  diss. 

empirical 

ch4 

1 

1 

1 

1 

1 

1 

c2h4 

0.278 

0.323 

0.311 

0.311 

— 

0.270 

(  +  16%) 

(  +  12%) 

(  +  12%) 

(-3%) 

C3Hs 

0.789 

0.810 

0.788 

0.788 

1.05 

0.749 

(+3%) 

(+0%) 

(+0%) 

(+33%) 

(-5%) 

h2 

0.0124 

0.0404 

0.00896 

0.00896 

0.00891 

0.0162 

(+226%) 

(-28%) 

(-28%) 

(-28%) 

(+31%) 

4.4  Blowout  predictions  and  experimental  results 

By  assuming  that  blowout  occurs  when  the  flame  base  reaches  some  minimum  value 
of  r],  r/i,,  each  of  the  theories  can  be  extended  to  predict  a  blowout  velocity  Uf,.  This 


assumption  is  equivalent  to  the  criterion  that  hb  be  proportional  to  the  axial  distance  at 
which  1 ,  /  =  where  Y^tab  =  Yst  for  the  mixing,  strain-rate  and  scalar  dissipation-rate 

theories,  and  =  1)  for  the  turbulent  flame  speed  and  empirical  theories. 

Using  Eq.  (4)  to  eliminate  h  from  Eq.  (29), 

jr  =  {c„A/c2)(-^)2TCtnfn(iib)/ f2{r}b).  (32) 

C  6  to 

Taking  ?//,  to  lie  the  same  for  all  fuels  results  in  the  blowout  prediction 


d* 

Ub 


/  1  stab  ,2 

^n,3\  -r-  )  7”c,n* 

i  o 


(33) 


Values  of  d*/l'i,  can  be  compared  with  Y'2tahTc ,n  for  each  theory.  These  are  shown  in 
Table  4.  where  the  values  have  again  been  normalized  by  those  of  CH.\. 

Only  the  authors  of  the  mixing  time  theory  and  the  strain-rate  theory  have 
actually  proposed  this  scheme  for  predicting  blowout.  Pitts,  in  his  empirical  correlation 
development,  began  by  assuming  that  blowout  occurs  when  the  stabilization  contour 
reaches  its  maximum  radius  rmax\  this  is  equivalent  to  choosing  a  particular  value  of  rjb. 
He  then,  however,  found  an  empirical  expression  for  the  local  velocity  at  blowout  which 
is  inconsistent  with  Eqs.  (25)  and  (26).  Effectively,  he  eliminated  } ;  from  Eq.  (26.1)  in 
his  discussion  of  blowout  which  then  changed  rc  c  to  c„/S2.  In  the  current  discussion, 
the  expression  of  Tr,t  consistent  with  Eq.  (26.1)  is  used. 

All  of  the  theories  make  reasonable  predictions  about  blowout  except  in  the  case  of 
hydrogen  where  the  turbulent  flame  speed  and  empirical  theories  are  far  off.  This  is  due 
to  these  theories  using  the  assumption  that  Ystab  =  1/  rather  than  Y„t  as  the  remaining 
theories  assume.  For  hydrogen.  Y)  is  very  different  from  Yat,  and  the  choice  will  affect 
predictions  greatly. 


Table  4 


Comparison  of  blowout  predictions  and  experimental  results 


Fuel 

d*/U  6 

turb.  f.s. 

mix.  time 

strain 

scalar  diss. 

empirical 

ch4 

1 

1 

1 

1 

1 

1 

c2h4 

0.416 

0.407 

0.358 

0.358 

— 

0.341 

(-2%) 

(-14%) 

(-14%) 

(-18%) 

c3h8 

1.10 

0.950 

0.861 

0.861 

1.15 

0.879 

(-14%) 

(-22%) 

(-22%) 

(+5%) 

(-20%) 

h2 

0.00382 

0.0529 

0.00464 

0.00464 

0.00460 

0.0213 

(  +  1300%) 

(+21%) 

(+21%) 

(+20%) 

(+460%) 

5.  Conclusions 

Five  theories  predicting  turbulent  jet  flame  liftoff,  each  utilizing  self-similar 
properties  of  isothermal  jets,  have  been  investigated.  The  theories  are:  a  turbulent  flame 
speed  theory,  a  cascade  mixing  time  theory,  a  maximum  strain-rate  theory,  a  limiting 
scalar  dissipation-rate  theory  and  an  empirical  theory.  Each  predicts  a  relation  between 
h  and  flow  variables  of  the  form 


h  Cn,\(^  stab /Yo)Tc,nt/ o  f n{T])i 


(29) 


with  different  expressions  for  Ystai,,  rc>n  and  fn(r]).  By  using  the  equation  for  the 
mean  concentration  field  and  by  expressing  the  similarity  functions  as  Gaussian,  rj  was 
eliminated  from  Eq.  (29)  to  arrive  at  the  relation 


b stab  h 

~yTT* 


V 

„  u  1  stab  \i 

cn,2\\  v  ) 


Uo 


o  \ 

c’n  d* ) 


En 


(30.2) 


(The  one-dimensional  cascade  mixing  time  theory  did  not  require  the  assumption  of 
Gaussian  functions  to  arrive  at  Eq.  (30.2).) 


The  five  theories  were  compared  for  their  ability  to  correctly  predict  the 
experimentally  observed  linear  dependence  of  h  on  UQ  by  examining  the  exponent  En. 
They  were  also  compared  for  their  ability  to  predict  the  relative  liftoff  heights  of  different 
fuels  by  examining  the  expression  for  Y3tabTc,n-  Finally,  the  theories  were  extended  to 
predict  blowout  and  compared  to  experimental  results. 

Because  the  theories  have  much  in  common  it  is  difficult  to  conclude  which  one  is 
the  best.  The  turbulent  flame  speed,  cascade  mixing  time  and  scalar  dissipation-rate 
theories  all  predict  highly  nonlin  ar  relations  between  h  and  Ua.  The  strain-rate  field 
in  turbulent  jets  has  not  been  measured  and  as  a  result  the  strain-rate  theory  cannot 
be  properly  evaluated.  The  empirical  correlation  has  no  physical  basis.  The  predictions 
of  relative  liftoff  heights  of  different  fuels  were  all  within  ±30%  for  all  of  the  theories 
and  all  of  the  fuels  except  the  turbulent  flame  speed’s  prediction  for  hydrogen.  Blowout 
predictions  of  all  the  theories  are  within  ±25%  except  for  the  turbulent  flame  speed  and 
empirical  theories’  predictions  for  hydrogen. 

The  similarities  in  the  theories  illustrate  their  limitations  in  understanding  lifted 
diffusion  flames.  To  make  significant  progress  in  understanding  the  liftoff  mechanism, 
more  detailed  physical  models  need  to  be  developed. 
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